Thermal rectification in bulk materials with asymmetric shape 
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We investigate thermal rectification in a bulk material with a pyramid shape to elucidate shape 
dependence of the thermal rectification, and find that rectifying coefficient R is 1.35 for this shape, 
which is smaller than R—1A3 for a rectangular shape. This result is fully duplicated by our numerical 
calculation based on Fourier's law. We also apply this calculation to a given shape, and show a 
possible way to increase R depending on the shape. 



A thermal rectifier is a device in which heat flows in 
a forward direction, while it hardly flows in the opposite 
direction. Owing to controllability of the heat current, 
the thermal rectifier can be applied to phononic devices 
as a diode is essential for electronic devices. Since the 
first report on thermal rectification at an interface be- 
tween Cu and CuO [1], this research field has been devel- 
oped by both experimental and theoretical studies [2] . In 
particular, several advanced studies concerning a discov- 
ery of thermal rectification in non-uniformly-mass-loaded 
carbon nanotube [3] and a prediction of thermal rectifi- 
cation in anharmonic ID lattice model/quantum thermal 
systems [4-7] have been reported in 2000s, which gives re- 
newal interests to this research field. After these studies, 
thermal transistor jsi], thermal memory [9|], and thermal 
logic gate [10^ are proposed as possible applications of 
the thermal rectifier. 

Recently, we have demonstrated thermal rectification 
of an oxide thermal rectifier as a bulk effect [Tl] accord- 
ing to a theory based on Fourier's law [l2|. This theory 
predicts thermal rectification in a sample made of two dif- 
ferent materials with different temperature dependences 
of thermal conductivities. According to Fourier's law. 
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FIG. 1: (Color online) (a) Schematic figure of the thermal 
rectifier (LSCO-top configuration). (l)-(4) represent chromel- 
constantan thermocouples, (5) resistive heater. S{x) indi- 
cates sectional area of the thermal rectifier, (b) Tempera- 
ture dependence of the thermal conductivity of LaCoOs and 
Lao.rSro.sCoOs polycrystalline samples. 
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heat flux J is proportional to thermal conductivity hi. 
Since is a. function of position x and absolute tem- 
perature T(x), the averaged thermal conductivity /^av of 
the sample may change depending on a direction of the 
temperature gradient under finite temperature difference 
AT (See Fig. 1 of Ref. [U] or Ref. [12]). Thus, the heat 
flux in a forward direction (Jj) can be different from 
that in the reverse direction (J^) through Eq. (1) lead- 
ing to > 1, where R represents the thermal rectifying 
coefficient defined by a ratio of the heat fluxs in the for- 
ward and reverse directions {R = jj^)- We have found 
R = 1A3 at 40 K with AT = 60 K in the thermal rectifier 
made of LaCoOs and Lao.rSro.aCoOs [11]. 
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There are several strategies to achieve larger R in bulk 
materials as follows. (1) To utilize an interface between 
two materials with different phonon bands. In this case, 
R is predicted to be up to 2000 [13]. (2) To provide mate- 
rials which exhibits structural phase transition accompa- 
nying a rapid change of at the transition temperature. 
(3) To form a thermal rectifier with a proper shape. It is 
theoretically predicted that carbon nanohorn and carbon 
nanoribon exhibit thermal rectification because of larger 
heat current from a larger section to the smaller section 
p!3-[l8j. This kind of the shape effect is also possible in 
bulk materials according to Eq. (1). However, experi- 
mental studies on shape-dependent thermal rectification 
have not been reported yet. In this letter, we will address 
the issue(3). 

Polycrystalline samples of LaCoOs (LCO) and 
LacrSrcaCoOa (LSCO) were prepared by a solid state 
reaction. Stoichiometric amounts of La203, SrCOa and 
C03O4 were thoroughly mixed, and calcined at 1273 K 
for 24h. The products were finely ground, filled up in la- 
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FIG. 2: (Color online) Monitored temperatures as a function 
of time at points (l)-(4) of (a) LCO-top and (b) LSCO-top 
configurations, respectively, (c) Temperature distributions in 
the thermal rectifier with LCO-top and LSCO-top configu- 
ration, respectively. Solid and dotted lines represent numer- 
ically calculated temperature distribution by Eq. (2). S{x) 
and l(x) are sectional area and the length of periphery of S{x), 
respectively. 



tex tubes with a small amount of water to increase sam- 
ple's porosity, and pressed at 40 MPa. Then, the pressed 
sample was calcined at 1273 K for 24h. X-ray diffraction 
pattern was measured using a standard diffractometer 
with Fe Ka radiation. Any impurity peaks were not de- 
tected. A mass density of the sample is evaluated to be 
50.7% of the ideal density for LCD, and 54.5% for LSCO, 
respectively. 

Our measurement system is depicted as a schematic 
figure in Fig.l (a). First, LCO and LSCO sample bars 
were cut to let them have a pyramid-like shape. The sec- 
tional area and length of the samples are shown in the 
right panel of Fig. 1(a). Then, the samples were bonded 
by silver paste with high thermal conductivity (Kyocera 
Chemical CT285), and annealed at 423 K for 0.5 h and 
at 493 K for 1.5 h to dry the paste. To make heat cur- 
rent /h into the rectifier, a resistive heater with r = 120 
Q was put on the top of the rectifier together with a 
copper plate at point (5). Another copper plate was at- 
tached on the bottom of the rectifier and the plate was 
glued to the cold head of closed cycle refrigerator by GE 
varnish (GE7031). Chromel-constantan thermocouples 
were attached at points (l)-(4) to monitor temperatures 
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FIG. 3: (Color online) Calculated and measured (a) heat cur- 
rents and (b) the rectifying coefficient R of both configura- 
tions. 



in the presence of the heat current. Detailed descrip- 
tion of the measurement system is given in the previ- 
ous reports [Tl], |l9| . Thermal conductivity was measured 
by a conventional steady state technique using chromel- 
constantan differential thermocouple. As shown in Fig. 
1(b), for LCO and LSCO is 5 times smaller than that of 
a single crystal [20], respectively, which is attributed to 
the high porosity of the samples. This low thermal con- 
ductivity is important to maintain the large temperature 
difference between points (1) and (4). 

Fig. 2(a) and 2(b) show the monitored temperatures as 
a function of time at points (l)-(4) in the forward direc- 
tion (LSCO-top configuration) and the reverse direction 
(LCO-top configuration), respectively. With increasing 
electric current /, the temperatures systematically in- 
crease. The heat current /h caused by the electric current 
through the relationship of /h = rP is controlled to keep 
70, 80, 90, and 100 K at point (1). The most striking 
feature is that the magnitude of the heat current in the 
forward direction (/hf ) is significantly larger than that in 
the reverse direction (/hr) showing the thermal rectifica- 
tion. The temperature difference between points (3) and 
(4) is completely explained by the thermal conductances 
of LCO and LSCO, which ensures that contact thermal 
resistance at the interface between LCO and LSCO is 
negligible. Note that contact thermal resistance between 
the rectifier and the cold head is also negligibly small. 

We evaluate this experimental results using Fourier's 
law. From Eq. (1), the temperature distribution of the 
thermal rectifier is derived as 



T{x) = T(0) 



4(0 



Jh(x) = Jh(0) -a r mmO' - T^}di, 
Jo 



(2) 



where /h, -S, <j, and Tc represent heat current, sectional 
area, Stefan-Boltzmann constant, length of periphery of 
the sectional area 5, and temperature of the cold head (40 
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FIG. 4: (Color online) (a) Two equivalent thermal rectifiers 
which show same calculation result of /hf(0), /hr(0) and R un- 
der same boundary condition, (b) Calculated as a function 
of X and y. The inset shows the enlarged calculated result in 
the range of 1 < x < 1.5, 1.425 <R< 1.465. 



K), respectively. The radiation loss is properly treated in 
the second line of Eq. (2). It should be emphasized that 
/h(0) is only parameter for calculating the temperature 
distribution under a boundary condition decided by the 
experiments. Here, as the boundary condition, temper- 
atures at point (1) {x = mm) and point (4) {x = 1.9 
mm) were set at 100 K and 40 K, respectively. As hz[x. 
T(x)], S{x)^ and /(x), experimental data shown in Fig. 1 
were used. As shown in Fig. 2 (c), calculated temperature 
distribution of the rectifier is in good agreement with the 
experimental data, which shows that the heat current in 
the pyramid-shaped rectifier is consistently explained by 
Fourier's law. 

Fig. 3 (a) shows the calculated and experimental 
heat currents /hf (0) (LSCO-top configuration) and /hr(0) 
(LCO-top configuration) as a function of temperature 
difference AT = T{x =0 mm)-T{x =1.9 mm). With 
increasing AT, /h monotonically increases, and experi- 
mental data are well reproduced by the calculation. As 
shown in Fig. 3 (b), the rectifying coefficient as a ratio 
of /hf (0) and /hf (0) is 1.35 at AT = 60 K, which is rather 
small compared with R = 1.43 of the rectangular one. 



To address how R depends on the thermal rectifier's 
shape, we perform numerical calculation of /? in a sam- 
ple made of two different homogeneous materials with " a 
given shape". Fig. 4 (a) shows two equivalent thermal 
rectifiers which show same calculation results of /hf, /hr, 
interface temperature (Tint) between LCO and LSCO, 
and R under same boundary condition. To convert from 
a pyramid shape (a given shape) to an equivalent rect- 
angular shape, a simple analytical relation is applied. 



Jo 



dx 
Six) 



(3) 



where L{L')^ and S{S') represent a length and a sectional 
area of a given shape (equivalent rectangular shape), re- 
spectively. 

Figure 4(b) shows the calculated as a function of a 
('^LSCO' '^LCO- sectional areas of 



and /3, where a 
LSCO and LCO with the equivalent rectangular shape). 
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lengths of LSCO and 



LCO with the equivalent rectangular shape). In the case 
of /3 = 1, takes the maximum of 1.46 at a = 1.2. Thus, 
R increases by 0.7% compared with a rectifier with a =1 
and (3 =1. Depending on R profile shifts on x axis, 
which is easily understood by thinking the thermal con- 
ductance of the thermal rectifier with rectangular shape. 
By applying this study to other materials, the thermal 
rectification may further enhance. We believe that this 
study will give a possible route to increase rectifying co- 
efficient. 

We have investigated shape-dependent thermal rec- 
tification from both experimental and theoretical ap- 
proaches. The cobalt oxide thermal rectifier with 
pyramid-like shape exhibits the rectifying coefficient 
of 1.35, which is rather small compared with that of 
rectangular-shaped thermal rectifier. The temperature 
distribution, heat current, and rectifying coefficient are 
consistently reproduced by Fourier's law. We have also 
shown that equivalent thermal rectifier with rectangular 
shape can be derived by a simple analytical relation, and 
found tiny enhanced rectifying coefficient at a = 1.2 and 
f3 =1. We believe that this study will give a possible 
route to increase rectifying coefficient. 

This study was partly supported by the Murata Sci- 
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